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Abstract : The importance of the study to the propagation of a Tsunami wave became from the 
complex phenomenon that represent and its natural disasters which represent a major risk for 
populations. To model this phenomena, we will consider the BBM-BBM 1 Boussinesq system [2] 
with a variable bottom in space and apply this system, first, using a mesh generated through a 
photo of the Mediterranean sea, then using a mesh generated through an imported xyz bathymetry 
for the East part of the Mediterranean sea. 

We choose here to use FreeFem++ [10] software which simplify the construction of the domain in 
particular, one of the advantage of FreeFem++ is that we can build a mesh through a photo and we 
can easily export bathymetric data in order to consider more realistic simulations. 

1 Introduction 

We consider here the numerical simulation of the BBM-BBM Boussinesq system in 2D over a 
variable bottom 

Vt + V.((D + V )V)-bV' (D 2 \/r] t ) = 0, 

(1) 

V £ + V77+^V|V| 2 -<LD 2 AV t = 0, 

where b and d are positive parameters (in the sequel, they are equal to 1/6). This system is a type 
of Boussinesq systems derived in [1, 5] as approximations to the three-dimensional Euler equations 
describing irrotational free surface flow of an ideal fluid Del 3 which is limited from below by the 
bottom —D(x,y) and from above by the free surface elevation r](x,y,t) (cf. Figure 1). 




Figure 1: The domain Q. 



1. Benjamin, Bona and Mahony (BBM) 
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The variables in (1) are non-dimensional and unsealed : X = (x, y) G ^ and £ > are pro- 
portional to position along the channel and time, respectively, 77 = /y(X, £) is proportional to the 

deviation of the free surface from its rest position, V = V(X,t) = ( J = ( u i v ) T = { u '-> v ) 



( d • 

is proportional to the horizontal velocity of the fluid at some height, V* = I ^ % J is the gradient, 
V-(*; , ) = d £C * +c^ # is the divergence and A # = d xx * + is the laplacian. 



Remark. In our study, we suppose that 77 = (9(a), where the amplitude a is the difference between 
the surface of water and the zero level, and the wave length A = in addition we limit our- 

selves to the case where rj + D > (there is no dry zone) since we are in a big deep water wave regime. 



The paper is organized as follows: in Section 2, we present the discretization of (1) in space an in 
time. Then, in Section 3, we present a method to build a mesh from a photo using FreeFem++ and 
in Section 4, another method to build a mesh through an imported bathymetry with FreeFem++ 
was detailed. Finally, in Section 5 we simulate the propagation of a wave, that looks like a Tsunami 
generated by an earthquake, in the Mediterranean sea over the BBM-BBM systems (1) with a flat 
bottom and with a variable one in space using the mesh generated from a photo of the Mediterranean 
and from the xyz bathymetry of a part of the Mediterranean sea. 



2 Discretization of the Boussinesq system 

We let n be a convex, plane domain, let Th denote a regular, quasiuniform triangulation of 
Q with triangles of maximum size h < 1, let Vh = {vh G C°((l);vh\T G Pi (T), VT G Th} denote a 
finite-dimensional subspace of H 1 ^) = {u G L 2 {0) s.t. |^ , |^ G L 2 {0)} where Pi is the set of 
polynomials of R of degrees < 1 and let (•; •) denote the L 2 inner product on Q. 
Consider the weak formulation of the system (1), find r]h,Uh,Vh G Vh such that \/<j)h G Vh we have : 

(Vhu <t>h) + b(p 2 Vr]ht] V0/^ + • (u h \ v h ) + rj hx u h + rj h u hx + rj hy u h + rf h u hy ] </> h ^ = 0; 
(Id + dV (D 2 ) • V) u h u <t>hj + d(p 2 Vu h t] ^<t>hj + (Vhx + u h u hx + v h v hx ] M = 0; 
(Id + dV (£> 2 ) • V) Uht; (j) h ^ + d(D 2 Vv ht ] V(j) h ^ + + + v/i^; <M = 0. 

(2) 

We note also that when using a mesh adaptation method integrated in FreeFem++, a P2-finite 
element space is used for better resolution. This method gives comparable results with Pi-finite 
element space without mesh adaptation technique however we have to more refine in the domain. 



Remark. The time discretization of the problem is realized using an explicit second order Runge- 
Kutta where the study of it convergence in the flat bottom case is detailed in [6] . 



3 Mesh generated through an imported a photo with FreeFemH — \- 

We present here a method to build a mesh from a photo using Photoshop® and a FreeFem++ 
script made by F. Hecht [9]. 

Through Google Earth®, we can download pictures of the earth for areas of interest. In this 
section, we choose the picture of the Mediterranean sea. The Mediterranean is being very large, in 
order to have better resolution, we take pictures of many parts (cf. Figure 2) that subsequently 
assembled using Photoshop® so as to obtain a complete picture of the Mediterranean (cf. Figure 
3). 

Using also Photoshop®, we can eliminate the dry areas that circumvent the Mediterranean (cf. 
Figure 4). 
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Figure 2: The pictures of two parts of the Mediterranean sea. 




Figure 3: The Mediterranean sea after assembly with Photoshop®. 



Then we convert the jpg photo to a pgm photo which can be read by FreeFem++ using in a 
terminal window : 

convert Medit_sea.jpg Medit _sea . pgm 



In order to generate the mesh of the Mediterranean sea domain (cf. Figure 5), we use a FreeFem++ 
script made by F. Hecht which can be downloaded from the following link 

http : //www . f reef em . org/f f 2a3/Stic-FF2A3-2010/edp/Chesapeake/Chesapeake-mesh . edp. 
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Figure 4: The Mediterranean sea. 




Figure 5: The mesh generated by FreeFem++ for the Mediterranean sea. 

Because of problem of non-smoothing boundary, we met some difficulties in the generation of 
the mesh in some regions as in the Figure 6 when using mesh adaptation technique. On the other 
side, without using the mesh adaptation technique, we no longer have this problem with the mesh 
generation since the mesh is generated at the beginning. 

Recently, and in order to smooth the boundary, F. Hecht has write another script and the 
result are better as shown in Figure 7. This script could be downloaded from the following link: 
http : / /lamf a . u-picardie . f r/ sadaka/FreeFem++/generate_a_mesh_using_a_photo . zip, and 
for more details see also [8]. 
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Figure 6: Mesh problem using the script "Chesapeake-mesh. edp" . 




Figure 7: The mesh around Cyprus using the script of F. Hecht (at left : the old one and at right : 
the new one). 

4 Mesh generated through an imported xyz bathymetry 
with FreeFemH — h 

In order to consider more realistic case i.e. the Mediterranean sea with a bathymetry, O. Pantz 
writes a script [11] with FreeFem++ that allows, through the xyz bathymetry of a chosen area of 
the earth ( downloaded from the NOAA 2 website, cf. Figure 8), to generate the mesh of the area 
where the amplitude is zero while using a level-set method to smooth the boundary (cf. Figure 9). 
In addition, we can through this script, on the one hand use the mesh adaptation method without 
having problem of generation of the mesh and in the other hand, we can have different label for 
each part of the boundary, which facilitates the use of different types of boundary condition. 

Using the triangulate function in FreeFem++, we start by reading the xyz downloaded data 
, then we separate the positive part from the negative one of the bathymetry by a curve which 
will be the boundary of our domain, we use the level-set method in order to smooth the boundary 
(cf. Figure 9). The details of this script could be downloaded from the following link: http: 
//lamf a.u-picardie . f r/ sadaka/FreeFem++/generate_a_mesh_using_xyz_bathymetry . zip 



2 . http : / / www . ngdc . noaa . gov/ mgg/ gdas/ gd_de signagr id . html 
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Figure 9: The mesh around Cyprus (left) and the mesh generated (right) using the script of O. 
Pantz. 
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5 Numerical simulations 



In this section, we will use FreeFem++ in order to simulate the propagation of a wave, that 
looks like a Tsunami generated by an earthquake, in the Mediterranean sea over the BBM-BBM 
systems (1) with a flat bottom —D(x,y) = —1 and with a variable one in space using the mesh 
generated from a photo of the Mediterranean sea (5) and from the xyz bathymetry of a part of the 
Mediterranean sea (7). 



5.1 Propagation of a Tsunami in the Mediterranean sea with a flat bot- 
tom. 

Inspiring from [4], we simulate here, the propagation of a wave, that looks like a Tsunami 
generated by an earthquake, in the Mediterranean sea over the BBM-BBM system with a flat 
bottom —D(x,y) = —1, where Dirichlet boundary conditions are taken on the hole boundary of 
our domain for 77, u and v and as initial datas : 



^0 = 0.1-^1- 1/(1 + 10 3 -exp 



-70 + (y - PY) - 0.2 • (x - PX) 
10 



-13 • 10 4 + 0.2 • (y - PY) + (x- PX) • 10 3 
10072) 

u h0 (x,y) = v h0 (x,y) = 0. 
Here PX = 2270, PY = 500, the step time is At = 0.1. 



In the Figure 10, we represent the maximum of the propagation of the solution at t = 1000sec, 
for the same datas above using Pi-finite element space and without mesh adaptation technique. 




Figure 10: The maximum of the propagation of the solution of a Tsunami in the Mediterranean 
sea for t = 1000 sec. Video available on : 

http : / /lamf a . u-picardie . f r/ sadaka/ f ilms/tsunami_mediterranee_MAX_Pl . avi. 



On the other side, we present in the Figures 11 and 12 the propagation of the same wave, using 
a mesh adaptation method and a P2-finite element space for better resolution. This method gives 
comparable results with Pi -finite element space without mesh adaptation technique, where in the 
last case, we no longer have a problem with the mesh generation, however we have to more refine, 
this makes the computation very heavy. 
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Figure 11: Propagation of the solution of a Tsunami in the Mediterranean sea for t = {0, 100, 400}sec. 
The Tsunami is generated near Crete. 

Video available on : http://lamfa.u-picardie.fr/sadaka/films/tsunami_mediterranee.avi. 
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Figure 12: Propagation of the solution and the mesh of a Tsunami in the Mediterranean sea for 
t = {0,100,400}sec. 

Video available on : http://lamfa.u-picardie.fr/sadaka/films/tsunami_mediterranee_ 
mesh . avi. 
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5.2 Propagation of a Tsunami in the Mediterranean sea with realistic 
bathymetry datas. 

In order to consider more realistic case i.e. the Mediterranean sea with a bathymetry, the 
bottom —D(x, y) of our study domain will be the same of the one downloaded , but since we are in 
a big deep water wave regime i.e. rj(x, y, t) + D(x, y) > 0, all the points of the bathymetry that are 
greater than or equal to —10m will be equal to —10m (cf. Figure 13). 

We remark also that to download the bathymetry data of the NOAA website, we must know the 
degree of Latitude and of Longitude of our domain. The spherical shape of the earth was taken into 
account, even if it does not play significant role because of the small spatial scale of the experiments. 
So, since the radius of the earth near the equator is R = 6378, 137/cm, we take into account that : 

1 degree of Latitude = tt ■ R/180 = 111, 263/cm, 

and 1 degree of Longitude = cos(l degree of Latitude • 7r/180) • tt • R/180. 

We limit our study for this example by making a translation of our domain considering that 1 
degree of Latitude is equal to 1 degree of Longitude is equal to 100km. We use an homogeneous 
Dirichlet boundary conditions for 77^, Uh and Vh for the shoreline boundary, and an homogeneous 
Neumann boundary conditions 77^, Uh and Vh at the open sea. We use also Pi -finite element space 
and we start from the following initial datas : 

/ /x-33.5- 100\ 2 /?/ - 33.8- 100 \ 2 \ 
% o = 0.01.exp(^ J To ) J> 

u h0 (x,y) = v h0 (x,y) = 0. 




Figure 13: The bottom of the bathymetry downloaded from the NOAA website, (min = —2.842km 
and max = — 0.01km). 



We present in Figure 14 the propagation of a wave generated near Cyprus in the Mediterranean 
sea thats looks like the one generated through an earthquake. We present in the Figure 15 the 
maximum of the propagation of the solution for t = 230sec. As we can see also in the same figure 
that a mesh adaptation technique is used too without having any difficulties in mesh generation. 
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Figure 14: Propagation of a Tsunami in the Mediterranean sea, generated near Cyprus for t = 75 
and 225sec. 

Video available on : http://lamfa.u-picardie.fr/sadaka/films/propagation_tsunami_ 
Chypre_Pl .mpeg. 
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Figure 15: Left : mesh adaptation around Cyprus after the script of O. Pantz. 
Video available on : http://lamfa.u-picardie.fr/sadaka/films/chypre_adaptmesh.avi 
Right : the maximum of the propagation of the solution of a Tsunami in the Mediterranean sea for 
t = 230sec. 

Video available on : http://lamfa.u-picardie.fr/sadaka/films/max_solution_tsunami_ 
Chypre_Pl . avi. 



6 Conclusion and Outlook 

We show in this paper, the utilities of FreeFem++ for the Boussinesq systemby building the 
domain, on the one hand through a photo taken from Google Earth® and on the other hand 
through an xyz bathymetry downloaded from the NO A A website. Concerning the simulation of a 
tsunami in the Mediterranean sea, the digital computing environment that we developed allows the 
integration of realistic data (bathymetry and geography) in a relatively simple framework. 

As we see, the mesh generation through a picture presents some difficulties, other then smoothing 
the boundary, we are limited to takes only homogeneous Dirichlet boundary condition on the 
hole boundary of our domain for 77, u and v. In the other hand, for the simulation of a Tsunami 
near Cyprus, we remark that an artificial reflection come from the open sea after t = 230sec 
which is due to the Neumann boundary conditions. In this case we must take some absorbent 
boundary condition which is an open problem. We see also, near the shoreline some oscillation 
because here, we are not in the big deep water wave regime and we must solve another equations 
such as Shallow Water equations [7] . Another approach is to use the sponge layer conditions as in [3] . 
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